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ABSTRACT 

The total 0.45-2.0 keV luminosity of M101 is 3.1 x 10 39 ergs s" 1 , of which 2.2 x 10 39 ergs s^ 
is due to diffuse emission. Of the diffuse emission, no more than 6% can be due to unresolved 
point sources such as X-ray binaries, and ~ 11% is due to dwarf stars The d^use emission 
traces the spiral arms and is roughly correlated with the Ha and FU emission, 
distribution closely follows the optical profile. The bulk of the diffuse emission is characterized 
by a two-thermal-component spectrum with kT = 0.20,0.75 keV, and the rati :os * off 1 the emission 
measures of the two components is roughly constant as a function of both radius and sur 
brightness. The softer component has a sufficiently large covering factor that the bulk of t e 
emission is likely extra-planar. We find no evidence of an extended axisymmetric X-ray halo, 
suggesting that any such halo has a strength much smaller than current predictions. 

Subject headings: X-rays-diffuse background — X-rays:ISM 


1. Introduction 

M101 (NGC 5457) is a nearby face-on spiral 
galaxy that lies in a direction of low Galactic ab- 
sorption. As such, it is an ideal target for the 
study of the hot component of the ISM in a “typi- 
cal” spiral galaxy, allowing high spatial resolution 
at relatively soft energies. 

1 Mailing address: Laboratory for High Energy Astro- 
physics, Code 662, NASA/GSFC, Greenbelt, MD 20771 

2 Universities Space Research Association 


1.1. Previous Results 

Previous studies of the diffuse emission from 
M101 have been executed with Einstein (Mc- 
Cammon k Sanders 1984; Trinchieri et al. 1990), 
ROSAT PSPC (Snowden k Pietsch 1995; Wang 
et al. 1999), ROSAT HRI (Wang et al. 1999), and 
ASCA (Wang et al. 1999). These studies are char- 
acterized by relatively poor spatial resolution that 
hindered point source detection and placed large 
uncertainties on the amount of diffuse emission. 

The PSPC study demonstrated the existence of 
diffuse j keV emission with a filling factor of ~ 1 
and a color temperature of 10 5 ' 8 K (kT = 0.054 
keV) The temperature and large filling factor 
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Table 1 

M101 (NGC 5457) Particulars 


Parameter 


Distance 

Position 

Inclination 
Line of Nodes 
N/f (foreground) 
Hubble Type 
Hi Mass 
Mass 

D25 


M101 

^ a ^ ue Source 


7.2 Mpc (2.09 kpc/arcmin) 
a = 14 /l 0I m 26.6- s 
6 = 54°35'19" (2000) 

27° ±1° 

35° ± 1° 

8.8 x 10 19 cm" 2 
Scd 

2.1 x 10 l ° M 0 
2.1 x 10 11 M 0 
1.0 x 10 n M 0 
23. '8 (49.8 kpc) 


Stetson et al. (1998) 


Kamphuis (1993) 

Zaritsky et al. (1990) c 
Hartmann & Burton (1997) 
Tully (1988) 

Kamphuis (1993) 

Kamphuis (1993) 

Comte et al. (1979) . 

Tully (1988) 


v c 

L fir 
SN Rate 


~ 170 km S" 1 

(3.4 ± 0.3) x 10 43 erg s~* ld 

2 century" 1 


Kamphuis (1993) 

Beck & Golla (1988) 
Matonick & Fesen (1997) 


Milky Way 

Value Source 


Sb-Sbc 

3.5-7.0X 10 9 M 0 
7 i 2.5 x 10 u M 0 
2A± l 0 3 x 10 u 
24.63 kpc 
26.71 kpc 
~ 220 km s“ 1 
3.0 x 10 43 erg s" 1 


Kerr (1993) 

Dame (1992) 

Binney & Merrifield (1998) 
Little & Tremaine (1987) 
Binney & Merrifield (1998) 
van der Kruit (1990) 

Sodroski et al. (1997) 


val ™* Zantsky e * al ( 199 0) values are representative of a number of measures made in the 1970s and 1980s, though the spread of 
values is somewhat larger than the uncertainties would indicate. p 

d Value summed over the entire disk. Note that the value in Fullmer & Lonsdale (1989) is for a point source at the location of M101. 


prompted the authors to suggest that the A keV 
emission was primarily in the halo; if the emission 
were from the disk and due to regions like the Lo- 
cal Hot Bubble (Snowden et al. 1998), with a typi- 
cal overburden of a few times 10 20 cm"* 2 , the filling 
fraction would be significantly more than unity, 
and the temperature would be much smaller than 
that observed for hot Galactic bubbles. Given 
a resolution of ~ l 7 , a primary concern was the 
contribution to the diffuse emission by unresolved 
point sources. Arguments were made that there 
was no known population of sources that could 
produce the amount of | keV emission seen, 

A PSPC study of the £ keV band (R45, 0.44- 
1.21 keV) is more comparable to the present Chan- 
dra study. Snowden & Pietsch (1995) determined 
an R45 band color temperature of 10 6 53 K (kT = 
0.29 keV), a flux within the inner 7.'5 of 5.8 x 10“ 13 
ergs cm" 2 s _1 , and a luminosity of 3.6 x 10 39 ergs 
s' 1 . Again there was a concern about the contri- 
butions from unresolved point sources, but there 
were insufficient source statistics to attempt a cor- 
rection. Contributions from sources such as M 
dwarfs could be reasonably well constrained and 
shown to be small. 


Wang et al. (1999) also came to the conclusion 
that no known population of sources could explain 
the entire emission seen in the central region of 
M101. They further attempted to constrain the 
temperature of the diffuse emission in the central 
5' with a joint fit to the ROSAT PSPC and A SC A 
GIS spectra. Given the angular resolution of the 
PSPC, and the poorer angular resolution of the 
GIS, they did not attempt to remove the point 
sources, but attempted to compensate by model- 
ing the aggregate spectrum of the point sources. 
Modeling the point sources with either thermal or 
power law spectra produced a diffuse component 
with kT = 0.186 keV (T=10 6335 ). 

1.2. Aim of the Present Work 

The angular resolution of Chandra not only al- 
lows point source detection and removal to un- 
precedented levels, but also allows the construc- 
tion of the point source luminosity function, the 
extrapolation of which places strict upper limits 
on the possible contribution of unresolved sources 
to the diffuse emission. The catalogue of point 
sources and the luminosity function were pre- 
sented in Pence et al. (2001), while discussion of 
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some individual objects is contained in Snowden 
et al. (2001) and Mukai et al. (2002). The purpose 
of the present work is to determine the distribu- 
tion, luminosity, and gross spectral characteristics 
of the soft diffuse X-ray emission from Ml 01. 

Due to the electronic noise feature, the Chandra 
energy interval is approximately E > 0.45 keV. 
In this energy interval the disk of M101 ranges 
from nearly transparent at all energies to, at worst, 
r ~ 2.7 at 0.45 keV. As M101 is nearly face on, we 
are “X-raying” the entire disk, detecting most of 
the gas hot enough to produce significant emission 
above 500 eV. The observed emission will contain 
contributions from HU regions, SNR, hot bubbles, 
and galactic chimneys, but the faintness of each in- 
dividual object, or the lack of contrast from neigh- 
boring emission regions, restrict the spectral anal- 
ysis to the bulk properties of the hot gas. 

2. Data 

The data analyzed here are derived from a sin- 
gle 98.2 ks Chandra exposure taken in AOl us- 
ing the ACIS-S3 back-side illuminated chip. The 
aimpoint was placed near the nucleus and the roll- 
angle adjusted so that the bulk of the j keV emis- 
sion detected by ROSAT, as well as the giant H II 
regions NGC 5461 and NGC 5462, would fall on 
the S3 chip. The configuration was also chosen to 
maximize the radial extent over which we could 
study the emission, as well as to leave space for 
“blank” sky from which to determine the back- 
ground. Other chips were on during the exposure 
but, due to their distance from the aim-point or 
poorer low-energy response, are of little use for 
this analysis. The most recent gain files were 
used; response and redistribution functions were 
extracted from the CIAO Caldb 2.12. 

The light-curve was constructed for the entire 
chip in the 2.0-7.0 keV band to allow the detec- 
tion and removal of background flares. Three mi- 
nor flares were detected and discarded, reducing 
the exposure to 92.3 ks. There are often low-level 
flares that are undetectable from the light-curve, 
but which produce a significant spectral signature 
in the 2.0-7.0 keV band. The contribution from 
these flares is discussed in §3.1.1. 


2.1. Point Source Removal 

Although the bulk of the diffuse emission emits 
at energies less than ~ 1500 eV, in order to pro- 
duce the cleanest spectrum of the diffuse emission, 
we removed sources detected to 3.5 <r in any of 
a series of energy bands covering the the entire 
Chandra energy range. Point sources were de- 
tected using the wavdetect routine in the 350-1300 
eV, 500-2000 eV, and 2000-8000 eV bands. (All 
of the sub-bands used in Pence et al. (2001) were 
processed in a similar way, but did not produce 
significant sources that were not found in three 
bands listed here.) Significance was determined 
using Poisson statistics; the source counts being 
determined from an elliptical 75% encircled energy 
region with the same aspect ratio and orientation 
as the PSF, and the background counts being de- 
termined from a concentric elliptical annulus ex- 
tending from 2 to 3 times the 95% encircled energy 
radius. A larger region was removed from around 
source P98 (Pence et al. 2001) due to its extreme 
brightness. 

The log N-log S relation in the 0.5-2.0 keV en- 
ergy interval is shown in Pence et al. (2001) and 
extends to luminosities ~ 10 36 erg s _1 . Extrapo- 
lating the log N-log S to lower luminosities is not 
trivial due to rapid changes in the source popula- 
tion with luminosity. In the Galaxy, the Low-Mass 
X-ray Binary (LMXB) population dominates for 
1036 erg s — 3 <; L ^ 10 33 erg s 3 , but tails off 
sharply below 10 36 erg s _1 . High-Mass X-ray Bi- 
naries (HMXB) with evolved companions have lu- 
minosities 10 35 erg s - 1 < L < 10 37 erg s" 1 , but 
only half the spatial density of persistent LMXBs. 
Below 10 35 erg s _1 , a number of other popula- 
tions exist, such as persistent Be X-ray binaries, 
whose density with respect to that of LMXBs is 
poorly understood. Assuming a source spectrum 
with T = 2 and directly extrapolating the Pence 
et al. (2001) log N-log S relation to luminosities 
10 35 and 10 34 erg s -1 produces 4.1% and 5.8% of 
the observed diffuse emission. 

3. Large-Scale Diffuse Emission 

3.1. Morphology & Photometry 

Figure 1 displays the chip S3 data after point 
source removal and smoothing by a 779 FWHM 
Gaussian. The lower parts of Figure 1 contain the 
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Fig. L- Top: The X-ray image in the 0.45-1.0 keV interval after point source removal, smoothed by a 
Gaussian with a HWHM of 3'.'94 (8 pixels). The color intervals are 1.6 x 10' 6 counts cm' 2 s' 1 arcmin' 2 
starting at the value of the background. White regions surrounded by high surface brightness regions denote 

^ en ‘° va so f rces Wlth lar S e p SFs. The irregular polygons are the regions used to determine 

the background. The rule is marked in arcminutes from the nucleus. Bottom Left: Optical image from the 
Palomar Sky Survey, with the position of the S3 chip marked. The ellipse marks the size and position of 

the arge-scale CO bar. Bottom Right: HI 21 cm image. The color scales runs from 10 2 ° cm' 2 (white) to 
2 x 10 cm . ' 
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optical and HI 21 cm images. It is immediately 
apparent that the X-ray emission is strongly cor- 
related with the spiral arms, and that a significant 
fraction of the X-ray emission is due to the giant 
HII regions NGC 5461 and NGC 5462. The nu- 
clear emission is, in fact, smaller than the emission 
in either of the giant HII regions. Many, but not 
all, of the knots of X-ray emission are correlated 
with HII regions (see Figure 2), and the stronger 
knots have associated “point sources” which may 
be, alternatively, compact diffuse emission regions 
or compact clusters of X-ray sources. The corre- 
lation of X-ray emission with the HI holes is as 
poor as the correlation between HI holes and Ha 
emission noted by Kamphuis (1993). (Since the 
lifetime over which an HI hole can be observed is 
~ 5 x 10 7 years, and the generating HII regions 
last ~ 2 x 10 7 years, Kamphuis argued that two- 
fifths of the holes should contain HII regions. In 
our sample, there are thirteen large HI holes, three 
are clearly filled with emission stronger than that 
of the surrounding disk, another three have X-ray 
emission comparable to the surrounding disk, and 
seven have less X-ray emission than the surround- 
ing disk. The HI holes that are filled by X-ray 
emission tend to have HII regions that are not 
near the peak of the X-ray emission, while X-ray 
empty HI holes tend to have X-ray emission and 
correlated HII regions on their edges.) 

We have compared the Chandra image with a 
FUV ( 152 lA) image from the Ultraviolet Imag- 
ing Telescope (UIT) (Stecher et al. 1997). As can 
be seen in Figure ??, the regions of diffuse X-ray 
emission are in general well correlated with the re- 
gions of diffuse UV emission. Most knots of X-ray 
emission are correlated with knots of U V emission, 
though there are some notable exceptions particu- 
larly the diffuse emission around P98, which may 
in fact be due to the wings of the PSF, and the 
arc-like emission region north-east of the nucleus. 
A number of the more intense FUV emission knots 
are not associated with peaks in the X-ray emis- 
sion. Some of these knots of FUV emission are on 
the edges of knots of X-ray emission, perhaps sug- 
gesting the propagation of star formation, the knot 
of X-ray emission being associated with SNR or 
an aggregation of young X-ray binaries, while the 
adjacent FUV emission region is associated with 
young high-mass stars. 


3.1.1. Determining the Background 

We had placed the nucleus near the aim-point 
of the S3 chip and specified a roll angle so that 
the bulk of the ROSAT £ keV emission would fill 
the S3 chip. The expectation was that the dif- 
fuse emission detected by Chandra , being signifi- 
cantly harder would, like the ROSAT f keV emis- 
sion, cover a significantly smaller area, and that 
the corners would be relatively free of contamina- 
tion by M101. Instead, we found that Chandra's 
greater sensitivity allows one to trace that emis- 
sion at least to the edge of the chip in almost every 
direction. 

In Figure 1 the three irregular polygons indi- 
cate regions with “low” surface brightness between 
strong spiral arms. The three regions contain sta- 
tistically identical spectra. The mean spectrum 
from these regions, after the removal of the instru- 
mental background (§3.2) is shown in Figure 3. 
There is significant flux in the 2-5 keV interval 
essentially all of which is due to the particle back- 
ground from unremoved low-level flares. By com- 
paring the background spectrum in the observa- 
tion of GK Per in periods with and without flares, 
one can show that the particle background during 
flares has a power law shape, without convolution 
with the response function. We have assumed that 
the spectral shape of low-level flares is similar to 
that of the stronger flares. To remove the residual 
low-level flares from the M101 observation, the 2- 
5 keV interval was fit with an unconvolved power 
law, and those parameters transferred to the fits 
in the 0.45-2.0 keV interval. 

The dashed line in Figure 3 is a model for 
the Galactic foreground/extragalactic background 
calculated for a 5° x5° box in the direction of M 101 
from the ROSAT All-Sky Survey (RASS) using 
the method of Kuntz & Snowden (2000). The 
extragalactic background (dotted line) is modeled 
by a power law of index 1.46 and a normalization 
of 10.5 keV cm" 2 s" 1 sr' 1 keV' 1 (Chen et al. 
1997). From the extragalactic background was re- 
moved the flux due to resoved point sources, as 
modeled by Mushotzky et al. (2000). The extra- 
galactic emission was then absorbed by the disk 
of M101 (~ 4.4 x 10 20 cm" 2 ). The Galactic fore- 
ground model contains contributions from the Lo- 
cal Hot Bubble (unabsorbed) and a two thermal- 
component Galactic halo absorbed by 9 x 10 
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fig. Left: I he diffuse X-ray emission in M101 compared to the HII regions from Hodge et al (1990) 

(crosses) the m holes listed by Kamphuis (1993) (ellipses), and the supernova remnants of Matonick & 

Fesen (1997) (boxes). Right: The UIT image (greyscale) compared to the X-ray image (contours). 



Fig. 3. Spectrum in the background regions. 
The best ROSA 7 model and the best fit Chandra 
spectrum are shown by the solid line. The dotted 
line is contribution by the unresolved background 
sources. 


cm -2 . Since the Raymond & Smith model tem- 
perature for the Local Hot Bubble is well con- 
strained by a series of ROSA T observations (Snow- 
den et al. 1998; Kuntz & Snowden 2000; Snowden 
et al. 2000), that parameter was fixed to T=10 611 
in the RASS fits, but all the other parameters were 
allowed to vary. The solid line is the best fit of 
this model to the Chandra data for the low sur- 
face brightness regions. 

It is not known whether the harder Galactic 
halo component is truly Galactic, or whether it 
may be in part extragalactic. Assuming it is ex- 
tragalactic and absorbed by the Ml 01 disk reduces 
the transmitted flux by only 4%. 

The relevant point is that the ROSATmodel for 
the Galactic foregrounds in the direction of M101 
produces a spectrum similar to that seen in be- 
tween the arms of M101, with nearly the same 
normalization. The best fit model to the Chandra 
data produces 77.8 x 10~ 6 counts s -1 arcmin -2 
in /?05ATband R45, while the ROSAT All-Sky 
Survey contains (81.8 ± 4.2) x 10 -6 counts s -1 
arcmin -2 for a square degree surrounding but not 
including M101. Thus, the low-surface-brightness 
interarm regions from which these spectra were 
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Fig. 4 — The 0.45-1.0 keV X-ray surface bright- 
ness as a function of radius. The solid histogram 
excludes the giant HII regions; the dashed his- 
togram includes the giant HII regions. The solid 
line is the fitted exponential. The dotted line is the 
optical profile from Okamura et al. (1976) scaled 
to the same central surface brightness for the ex- 
ponential disk. 

derived appear to have very little contamination 
by emission from M101 itself. We have used this 
spectrum as the Galactic foreground spectrum. 

3.1.2. Surface Brightness and Total Flux 

The radial profile of the 0.45-1.0 keV emission 
is shown in Figure 4. As might be expected, the 
two giant HII regions produce large perturbations 
to the radial distribution. Without the contribu- 
tion of the giant HII regions, the X-ray exponential 
scale length is 5.24 ± 0.49 kpc, similar to the scale 
lengths in both the V band, Rv = 5.19 kpc, and 
the B band, R B = 5.40 kpc (Okamura et al. 1976). 

Excluding the giant HII regions, the total flux 
within a radius of 10 kpc is 1.5 x 10 counts cm 
s _1 in 0.45-1.0 keV. The giant HII regions provide 
an additional 0.24 x 10 -4 counts cm -2 s -1 in 0.45- 
1.0 keV. Since the flux is dominated by the low- 
surface-brightness regions we used the spectrum of 
the low surface brightness gas to find a conversion 
to physical units of 1.60 x 10“ 9 ergs count -1 , and 
thus a total 0.45-1.0 keV luminosity of 1.77 x 10 39 
ergs s -1 , or a 0.5-2.0 keV flux of 2.05 x 10 39 ergs 
s -1 . Within a radius of 10 kpc, the total luminos- 


ity of point sources (other than P98) is 0.82 x 10 39 
ergs s -1 . Extrapolating the exponential relation 
shown in Figure 4, assuming a T = 2.0 spectrum 
for the point sources, and adding the flux due to 
the nuclear region and the giant HII regions, we 
find a luminosity of 3.86 x 10 39 ergs s -1 in the 
central 7!5 in the 0.45-1.0 keV interval. This value 
compares well to the 3.6 x 10 39 ergs s -1 found by 
ROSAT in the same region in the 0.44-1.21 keV 
interval (Snowden & Pietsch 1995), further imply- 
ing that about 45% of the “diffuse” ROSAT flux 
was due to unresolved point sources. 

3.1.3 . Color Gradients 

A brief glance at the spectra in following sec- 
tions reveals that the emission in the 0.45-2.0 ke\ 
interval is dominated by OVII (0.560 keV), OVIII 
(0.650 keV), and a complex of iron lines in the 
.7-1.0 keV interval. For a first-order understand- 
ing of the diffuse emission we considered the (.70- 
1.0)/(.45-.70) hardness ratio. Because the count 
rate is low, and thus the counting statistics poor, 
we present the hardness ratio as functions of ra- 
dius and of surface brightness (Figure 5) rather 
than as a map. The hardness ratio is constant 
at low surface brightnesses (below ~ 10~ 5 counts 
cm" 2 arcmin" 2 ). In the surface brightness in- 
terval (1.3- 1.8) x 10" 5 counts cm" 2 s" 1 arcmin" 2 
the hardness ratio is much higher. However, the 
bulk of this emission comes from a limited number 
of compact, bright emission regions, some of which 
are identified as HII regions (Hodge et al. 1990). 
Above 2 x 10" 5 counts cm" 2 s _1 arcmin" 2 , the 
emission comes from the nucleus and the giant HII 
regions. There is a slight and statistically insignif- 
icant increase of hardness with surface brightness 
for the nucleus, and a greater increase (though still 
statistically insignificant) for the giant HII regions. 
There is no statistically significant trend of hard- 
ness with radius. 

3.2. Spectroscopy 

The background spectrum is a combination of 
the Galactic foregrounds, the extragalactic back- 
ground, and the instrumentation particle back- 
ground. The particle background spectrum was 
derived from the u dark moon data and scaled 
to this observation /region using the particle back- 
ground rate in the 10.5-14.5 keV region, where 
Chandra has an insignificant response to cosmic 
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Fig. 5. — Top: The hardness ratio as a function 
of radius. The solid line excludes the giant HII 
regions; the dashed line includes the giant HII re- 
gions. Bottom: The hardness ratio as a function 
of surface brightness. The lower histogram is more 
broadly binned, and is offset in color by -1. The 
upper histogram above 1.3 x 10~ 6 counts cm -2 
s 1 arcmin 2 is dominated by a small number of 
bright knots. 

X-rays. The Galactic foreground/extragalactic 
background spectrum was derived from the back- 
ground regions noted in §3.1.1. The ratio of the 
response in the region covered by MIOl’s diffuse 
emission to the response in the background region 
is energy dependent, and can vary from one to 
five percent over a hundred eV. Since the par- 
ticle background is patterned (Markevitch 2001) 



Fig. 6. — A raw spectrum of the diffuse emission. 
The lower spectrum is the background. Features 
at 1.2 and 1.38 keV (see §3.2) are apparent in the 
raw spectrum; the 1.2 keV feature may have been 
augmented by the background subtraction. 

the particle background spectrum will not scale 
with area. Therefore, we removed the appropriate 
particle background spectrum from the spectrum 
extracted from the background regions, scaled the 
result by the ratios of the reponses and areas, and 
added the correct particle background spectrum 
for the diffuse emission, propagating all uncertain- 
ties. Figure 6 shows a raw spectrum of the large- 
scale diffuse emission and its background, as com- 
puted by the method outlined above. 

Since the exposure was made on day 247 after 
launch, the accumulating absorption thought to 
be due to organics condensing on either the chips 
or the filters is negligible. 

We have used solar metallicities in the spec- 
tral fitting. The abundance gradient in M101 has 
been studied extensively using gas phase abun- 
dances in HII regions (Garnett et al. 1999; Garnett 
& Kennicut 1994; Garnett 1989; Torres-Piembert 
et al. 1989; Evans 1986), yielding either abundance 
trends (He, N, 0, S) or mean abundances (C, Ne, 
Ar). For the region covered by the Chandra obser- 
vation, the abundances are close to solar for these 
elements. The nucleosynthetic studies of Timmes 
et al. (1995) and Alibes et al. (2001), which con- 
tain compendia of Galactic abundance data in the 
form of X/Fe as a function of Fe/H, can be used 
to link the observed species to those of interest for 
X-ray spectroscopy, e.g, Ca, Si, Fe, and Ni. 

We have extracted spectra for three values of 
surface brightness, A: 1.6-4.8, B: 4. 8-8.0, and C: 
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Fig. 7 — Top: Spectrum A, the spectrum for re- 
gions with surface brightness 1.6-4. 8 x 10 6 counts 
cm - 2 s — 1 arcmin -2 . The solid line is the fit as- 
suming no internal absorption, the dotted line is 
the fit assuming half of the emission is behind the 
full absorbing column of the disk. The inset shows 
the region from which the spectrum was extracted. 
Middle: Spectrum B, the spectrum for regions 
with surface brightness 4. 8-8.0 x 10 6 counts cm 
s -1 arcmin -2 . BottoTn: Spectrum C, the spectrum 
for regions with surface brightness 8.0-11.2 xlO -6 
counts cm -2 s -1 arcmin -2 . 


8.0-11.2 xlO -6 counts cm" 2 s -1 arcmin -2 (where 
the values are above the nominal background). 
These spectra are shown in Figure 7. 

Spectra A and B are well fit with two thermal 
components, using MeKaL models (Liedahl et al. 
1995; Mewe et al. 1986, 1985). That the spectra 
are fit by two thermal components is not surpris- 
ing. In the 0.45-1.0 keV interval, there are three 
major features to be fit: the OVII 0.56 keV line, 
the OVIII 0.65 keV line, and the 0.7-1.0 keV com- 
plex of Fe lines. The Fe complex dominates the fit 
of the harder component, while the OVII/OVIII 
ratio dominates the fit of the softer. Our simplest 
model does not include the effect of absorption in- 
ternal to M101. A trial fit where the absorption 
for each component was allowed to vary indepen- 
dently found no absorption for the soft component. 
The absorption for the hard component was, un- 
derstandably, poorly determined, with the uncer- 
tainty being nearly as large as the fitted value, 
~ 3 x 10 21 i cm -2 . 

After fitting the two thermal components, there 
remained two further features at 1.2 and 1.38 keV, 
as well as a small amount of flux below 1.65 keV. 
Above 1.65 keV, the Si Kq line from the particle 
background, though well removed, substantially 
increases the uncertainty. Above 2.0 keV there is 
no significant flux. Examination of Figure 6 shows 
that, while the 1.2 keV feature may be augmented 
by poor background subtraction, neither feature is 
produced by the background subtraction process. 
Non-equilibrium models were no more adequate at 
fitting these lines than MeKaL models. Due to the 
small number of counts in this energy interval, the 
spatial distribution of this contribution can be de- 
termined only from very heavily smoothed images 
(FWHM~ 32"). The distribution of the 1.1-1.65 
keV emission is poorly correlated with the 0.45- 
1.0 keV emission, which contains the bulk of the 
galactic emission. Although 1.1-1.65 keV emission 
is typically found where the 0.45-1.0 keV emission 
is strong, it is also found where there is little 0.45- 
1.0 keV emission. The 2.5-5.0 keV emission should 
contain no galactic emission, and should be com- 
posed primarily of the particle background; the 
1.1-1.65 keV emission is not correlated with the 
2.5-5.0 keV either. To account for these features 
in fitting, they were modeled as unresolved lines; 
Gaussians with only instrumental broadening. 

In order to estimate the effects of internal ab- 





Fig. 8. — The total absorbing column. The con- 
tours are logarithmic from 10 2 ° to 3.2 x 10 21 cm' 2 
in steps of 10° 25 . The heavy lines show the loca- 
tion of the Chandra image and the area covered by 
the CO map. With the exception of the giant Hll 
regions, the bulk of the region containing X-ray 
emission is covered by the CO map. 

sorption on the spectral parameters, we fit the 
spectra both assuming no internal absorption and 
assuming that half of the emission is behind the 
disk absorption. The fit parameters are listed in 
Table 2. The map of total absorbing column was 
constructed from the 30" resolution Hi map of 
Kamphuis (1993) and a 55" resolution CO map 
kindly provided by T. Heifer. A map of the total 
absorbing column density is shown in Figure 8. 
The CO map does not cover the entire area ob- 
served by Chandra , but does cover the area con- 
taining the bulk of the emission. A model assum- 
ing that half of the emission is behind the ab- 
sorption due to the M101 disk typically reduces 
the temperature of the soft component by ~ .005 
keV, increases the total soft emission measure by 
~ 87%, reduces the temperature of the hard com- 
ponent by ~ .009 keV, and increases the total hard 
emission measure by ~ 60%. 

The temperature of the soft component is sim- 
ilar to that found by Wang et al. (1999) (0.186 
keV) from a joint ROSAT PSPC - A SC A GIS fit. 
Likely as a result of a factor of ~ 3 better spec- 
tral resolution, we do not see evidence for their 


hotter component (kT >5.1 keV), and fixing our 
hard component to kT = 5.1 keV produces a sig- 
nificantly poorer fit (A* 2 = 28). 

3.2.1. Joint Fits with ROSAT 

Given that ROSAT has a much stronger soft 
response, and softer components could contribute 
to the OVII emission. These softer component, not 
visible to Chandra could change the temperatures 
of the components fitted to the Chandra spectrum. 
We attempted to determine whether the Chandra 
fit was consistent with the ROSAT data. 

From the ROSAT data we extracted the region 
covered by the Chandra chip S3, but not includ- 
ing the giant Hll regions or a small region around 
the nucleus. The background was derived from an 
annulus between 25' and 50', beyond the observed 
M101 emission. The raw source spectrum con- 
tained 3769 counts, of which 2538 were due to the 
particle background and Galactic foreground. The 
comparable Chandra spectrum was created in two 
parts. First, we extracted and fit the spectrum for 
all of the regions containing point sources, exclud- 
ing source P98 (since P98 is brighter than all of the 
other point sources combined, and was clearly not 
visible during the ROSAT exposure). The point 
source contribution was reasonably well modeled 
by a powerlaw with T — 1.6 and N H = 4 x 10 20 
cm' 2 . Second, we extracted and fit the diffuse 
emission from the remainder of the chip, exclud- 
ing the point sources, the giant Hll regions, and 
a small region around the nucleus. For both the 
Chandra and ROSAT spectra, the contribution of 
the unresolved extragalactic background was cal- 
culated and removed. 

The combined model is shown in Figure 9. The 
normalization of the Chandra model above the 
carbon edge (primarily the contribution from the 
diffuse emission) appears to be slightly too high, 
but within the uncertainties of the individual data 
points. Below the carbon edge ( ROSAT bands 
R1 and R2, the keV” emission), the emission 
is split nearly evenly between point sources and 
the diffuse component. There remains significant 
emission below the carbon edge that is not ac- 
counted for by the Chandra model. In joint fits 
of the ROSAT and Chandra data, the tempera- 
ture of the soft Chandra component remained well 
constrained, but neither the very soft component 
seen only by ROSAT nor the harder Chandra com- 
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Table 2 

Spectral Fit Parameters 


Spectrum 

Counts a 

R 

kpc 

^nucleon 

10 20 cm -2 

A 

12907/5582 

5.64 

0.0 

0.0, 3. 2 

B 

4459/1231 

4.80 

0.0 

0.0,14.2 

C 

1442/315 

4.76 

0.0 

0.0,15.4 

A+B 

17366/6813 

5.49 

0.0 

0.0,13.1 

ROSAT c 

23872/16640 

7.03 

0.0 

0.0,10.6 

Ext.Nuc. 

794/276 

1.12 

0.0 

0.0,23.3 


T S 

EM S 

T H 

kT 

x 10" 2 cm" 6 pc 

kT 

.210 ± 0.015 

1.75 ± .17 

.741 ± .048 

.203 ± .013 

1.31 ± .11 

.730 ± .080 

.195 ± 0.010 

3.72 ± .28 

.754 ± .066 

.190 ± .009 

2.87 ± .20 

.752 ± .073 

.197 ± .021 

5.66 ± .67 

.607 ± .117 

.194 ± .018 

4.58 ± .63 

.599 ± 101 

.204 ± .010 

2.11 ± .14 

.746 ± .060 

.198 ± .012 

1.59 ± .11 

.732 ± 064 

.215 ± .028 

.745 ± .135 

.695 ± .210 

.210 ± .025 

.540 ± .082 

.691 ± 237 

.198 ± .026 

3.48 ±.54 

.706 ± 126 

.196 ± .019 

3.06 ± .48 

.711 ± .153 


EMh 

xl0~ 3 cm -6 pc 


5.94 ± 1.06 
3.46 ± .54 

13.3 ± .2 
7.77 ± .97 
17.6 ± .6 

9.95 ± .47 
7.56 ± .91 
4.23 ± .54 
1.35 ± 1.14 
.704 ± .576 
8.53 ±3.56 
4.74 ±2.84 


X* 


0.71 

0.71 

1.56 

1.55 

1.36 

1.37 
0.85 
0.86 
.55 
.55 
1.07 
1.10 


i/ 


97 

97 

85 

85 

43 

43 

97 

97 

99 

99 

32 

32 


-The first number is the number of raw counts in 0.45-2.0 keV, the second number is the calculated number of part.cle background 
counts in the same energy interval. 


b The column density of hydrogen nucleons was fixed. All absorption is exclusive of the 9 x 10 
c Region used for comparison with ROSAT. 


19 cm' 2 due to the Milky Way. 


ponent were well constrained. Assuming that the 
fit to the Chandra spectra was a good model of 
reality, we fitted the very soft component in the 
ROSAT data and found the 90% confidence in- 
terval to be kT < 0.053 ( logT = 5.79), and the 
68% confidence interval for the normalization to 
include zero. 

Given the small number of counts, the difficulty 
of background subtraction below the carbon edge, 
and the uncertainty of the ROSAT calibration at 
the extreme low energy boundary, one can confi- 
dently state only that there appears to be some 
excess emission below 0.35 keV . One can not com- 
pletely eliminate the possibility that the softer 
component seen by Chaudvo, does, in fact, con- 
tain contributions by even softer components. The 
similarity of the upper limit of the very soft com- 
ponent to the temperature derived by Snowden & 
Pietsch (1995) from the ROSAT R2/R1 band ra- 
tio is surprising, and likely coincidental, given that 
such a large portion of the \ keV emission is due 
to point sources (§3.1.2). 


4. Small-Scale Diffuse Emission 

4.1. Nucleus/Bulge 

M101 has a large-scale bar in CO (J=l-»0) 
emission (Kenney et al. 1991), which is marked 
in Figure 1. It has a position angle of 102° ± 4°, 
has a length of l!5 ± 0.'5, and appears to “an- 
chor” the main optically defined inner arms. At 
higher resolution (2 ;/ ), Ml 01 is a classic twin- 
peak” source, with strong CO emission from two 
compact regions flanking the nucleus where the 
flow of molecular gas down the bar crosses the ILR 
Kenney et al. (1992). The two peaks are separated 
by - 375 and have a position angle of roughly - 
10°. In the CO images there are filaments stretch- 
ing from the twin peaks outwards, roughly along 
the bar 3 . There is no CO concentration associated 
with the nucleus. 

A similar morphology is seen in Ha (Moody 
et al. 1995). The (presumed) nucleus is bright in 

3 Note that the position angle of the bar described in (Ken- 
ney et al. 1991) is rotated by ~ 25° from the bar described 
in Kenney et al. (1992). 
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Fig. 9. — The ROSAT spectrum of the region cov- 
ered by the Chandra S3 chip. The smooth solid 
line is the contribution by the unresolved extra- 
galactic background. The dotted line is the contri- 
bution by the point sources detected in the Chan- 
dra image. The dashed line is the sum of the Chan- 
dra model of the diffuse emission and the previ- 
ous two components. The stepped solid line is the 
previous spectrum with the addition of the fitted 
“ROSAT' component. The residuals plotted are 
from the fit containing the 11 ROSAT’ component. 
The grey residuals are from the Chandra model 
alone. 

Ha, and is flanked by two peaks of emission 2" 4 ± 
(XT5 from the nucleus (further out than the CO 
peaks) at a position angle of +10°. A strong Ha 
filament coincides with the southern CO filament, 
and some Ha features coincide with an extension 
of the northern CO filament. A WFPC I image 
(A5479) on PC chip 6 shows a bright point-like 
source at the position of the southern Ha peak, 
but no northern counterpart (Moody et al. 1995). 

Figure 10 shows the nuclear region of our Chan- 
dra observation smoothed by T/968 (4 ACIS pix- 
els) HWHM gaussian. Pence et al. (2001) used 
optical SNR to tie the X-ray coordinate system 
to the optical coordinate system of Matonick & 
Fesen (1997), and found that the Chandra coordi- 
nates need be corrected by an offset of 0'.'47±0'.'20, 
less than the size of an ACIS pixel. Boxes mark 
X-ray point sources; the position uncertainties are 
significantly smaller than the marked boxes. The 
separation of the two nuclear sources is S'/lb. 

The large cross is the position and uncertainty 



Fig. 10.— The nuclear region. The X-ray coun- 
tours are [1, 1.5, 2.0, 2.5, 3, 4, 5, 6, 7, 8]xl0“ 5 
counts cm -2 s _1 arcmin -2 . The image has been 
smoothed by a 3'.'97 FWHM Gaussian. The solid 
straight line is the line of nodes, the northwest cor- 
ner is towards the observer. The dashed straight 
line is the orientation of the molecular bar. The 
cross shows the location and uncertainty of the Is- 
rael et al. (1975) position of the radio continuum 
nucleus. The upper x is the position of the optical 
nucleus from a WFPC I image; the lower x is the 
position of the southern optical/Ha source. The 
small boxes are the X-ray point sources. The non- 
systematic uncertainties of the positions of the op- 
tical and X-ray sources is much smaller than the 
plot symbols. 

of the continuum nuclear source (Israel et al. 
1975), and the dynamical center determined from 
the CO (Kenney et al. 1991). The upper “x” is 
the optical nucleus and the lower “x” is the south- 
ern optical/Ha source. The non-systematic un- 
certainty for the optical sources is also very small 
on this figure. Both optical sources, the northern 
X-ray source, and the peak of the diffuse X-ray 
emission are within the positional uncertainty for 
the continuum nucleus. The long tails on the op- 
tical symbols show the shift necessary to make the 
optical nucleus coincide with the northern X-ray 
source; if the northern X-ray source is the opti- 
cal nucleus, the southern X-ray source does not 
coincide with the southern Ha source. 
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Since both' the northern and southern nuclear 
X-ray sources have spectra consistent with that 
of rather ordinary X-ray binaries, their locations 
may be coincidental. This region is dynamically 
interesting, and might be expected to shed light on 
the complex occurences at the center of our own 
Galaxy. If we take the HST and Chandra astrom- 
etry at face value, the bulk of the diffuse X-ray 
emission is centered near the southern optical/ Ha 
source, which, given its proximity to the CO peak, 
might be expected to be a region of intense star 
formation. Given the extended nature of the emis- 
sion and the position of the line of nodes, the lack 
of a northern counterpart is curious. Conversely, 
a ~ 175 shift, somewhat larger than the expected 
HST /Chandra r.m.s. offset (< 172), could make 
the optical nucleus coincident with the northern 
X-ray source. 

Our understanding of the nucleus is hampered 
by the lack of a good independent measure of the 
offset in the optical and X-ray coordinate systems. 
This measure will come with the execution of a 
survey of M101 with the HST ACS in Novem- 
ber 2002, which will attempt to identify the X-ray 
sources over the entire ACIS-S3 field with optical 
counterparts. 

4-1.1. Spectroscopy 

We extracted spectra from the nuclear region 
at two different surface brightness levels. From 
the compact nucleus itself we extracted^ the re- 
gion with surface brightness > 1.7 x 10 3 counts 
cm -2 s -1 arcmin -2 shown by the higher thick con- 
tour in Figure 11. The mean column density for 
this region is 2.8 x 10 21 cm- 2 , the bulk of which 
(2.6 x 10 21 cm -2 ) is due to molecular gas. Al- 
though the nucleus itself is quite bright, the to- 
tal number of counts is still small (239 counts in 
0.11 arcmin 2 ). The spectrum is reasonably fit by 
an absorbed Bremsstrahlung, N H = 1-8 x 10 21 
cm' 2 , kT = 0.262 keV, x* = 0.58 for v = 10 
(Figure lib). Other spectral forms (power law, 
thermal plasmas, etc.) produced worse fits. 

We extracted the spectrum of the extended 
emission around the nucleus, the region with sur- 
face brightness 0.725 x 10 5 counts cm 2 s 
arcmin -2 <c S <C 1.7 x 10 3 counts cm s 
arcmin -2 ; between the two thick contours in Fig- 
ure 11. The mean column density for this region 
is 2.3 x 10 21 cm -2 , the bulk of which (2.2 x 10 21 



Fig. 11.— Top: The extended nuclear region. 
The smoothing is the same as in Figure 1 (twice 
as much as in figure 10). The contours are 
from 0.32 x 10 -5 counts cm -2 s -1 arcmin 2 to 
3.2 x 10 -5 counts cm -2 s -1 arcmin -2 in steps of 
0.32 x 10 -5 counts cm -2 s -1 arcmin 2 . The dark 
contours at 0.725 and 1.7x10 5 counts cm s 
arcmin -2 delimit the regions from which spectra 
were extracted. Middle: The spectrum for the 
region with surface brightness greater than 1.7 x 
10 -5 counts cm -2 s -1 arcmin -2 . Bottom: The 
18 spectrum for the region with 7.25 x 10 6 <surface 
brightness< 1.7x10 5 counts cm 2 s arcmin 





cm -2 ) is due to molecular gas. This spectrum 
(Figure 11c) can be fit with a two thermal com- 
ponents (see ext.nuc. line in Table 2). Although 
we allowed the absorption to vary independently 
for the components, the fitted absorption was con- 
sistent with no absorption for either component. 
Including a Breinsstrahlung component as a con- 
taminant from the nucleus changes the thermal 
temperature and reduces the normalizations, but 
only slightly. 

Thus, as the nuclear emission appears to be ab- 
sorbed by about half of the total column density, 
it seems reasonable that the nucleus itself is em- 
bedded in the molecular ISM of the bar. The 
extended nuclear emission appears to be mostly 
unabsorbed, and is likely extraplanar. Assuming 
that half of the extended nuclear emission is be- 
hind the column density measured in Hi and CO 
does not significantly change the goodness of the 
fits or the fitted temperatures. 

4.2. Giant HII Regions 

Three giant HII regions fall on the S3 chip (see 
Figure 1). Of these, NGC 5453 was not detected 
by either R0SA7PSPC (Williams & Chu 1995) or 
HRI (Wang et al. 1999) surveys. In this Chandra 
image there is a very faint enhancement partially 
encircling the listed position. NGC 5462 fails on 
the part of the chip with the worst PSF, and the 
PSF for NGC 5461 is only slightly better. 

Both NGC 5461 and NGC 5462 were studied by 
Williams & Chu (1995) who found them resolved 
by the PSPC, but due to contamination by the 
soft extended disk emission they could not rule out 
the presence of a single large X-ray source produc- 
ing the X-ray emission from each region. Despite 
the relative softness of HII regions, both of these 
objects show extended emission in the deep HRI 
survey by Wang et al. (1999) In the Chandra im- 
age, both NGC 5461 and NGC 5462 show diffuse 
extended emission as well as bright point sources. 

NGC 5462: The spectrum was extracted from 
a region with a surface brightness greater than 
3.3 x 10 6 counts cm"" 2 s -1 arcmin"" 2 . The dif- 
fuse emission is reasonably fit (xl = 0.97) by an 
absorbed MeKaL model. A somewhat better fit 
( Xu = 1*2) can be found by assuming the sum of an 
unabsorbed thermal component and an absorbed 
thermal component, the same type of model that 
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Fig. 12. — Top: Spectrum of NGC 5462. The 
model spectrum is the sum of an unabsorbed 
MeKaL thermal component and an absorbed 
MeKaL component. The inset image contours are 
2.4, 4.0, 6.3, 10., 16., 24. xlO -6 counts cm -2 s -1 
arcmin -2 . The black regions are either off-chip or 
removed for point sources. The scale bar is 30" 
long. The box is the location of SN 1951H. Bot- 
tom: Spectrum of NGC 5461. The model spec- 
trum is an absorbed MeKaL thermal model. 
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fits the disk 'emission. (See Table 3.) In this 
case, the temperature of the unabsorbed compo- 
nent is well-constrained and similar to the soft 
disk component; the absorbed component is not 
well-constrained, but is clearly hotter. A fit con- 
taining an unabsorbed MeKaL component and an 
absorbed power law is significantly worse, imply- 
ing that the non-thermal component is small. Ex- 
tracting a spectrum using a higher surface bright- 
ness cut did not significantly change the fits. The 
supernova remnant SN 1951H is offset from the 
X-ray center, but is within the extended emission. 
Given that a giant HII region is likely to be formed 
of many different X-ray emitting components, ex- 
panding hot stellar wind bubbles, SNR, etc., and 
given the small number of counts in this spectrum, 
the poorness of the fit is not surprising. 

NGC 5461: The diffuse emission is reasonably 
fit {xl = 1.58) by an absorbed MeKaL model. As 
with NGC 5462, a somewhat better fit (xl = 1-2) 
can be found by assuming the sum of an unab- 
sorbed thermal component and an absorbed ther- 
mal component, but the fit parameters are poorly 
constrained. 

5. Discussion 

Observations such as those discussed in this pa- 
per stand at the curious intersection of disciplines: 
between the theory of galaxy formation and the 
classical study of the ISM, between processes dom- 
inated by the depth and shape of the dark matter 
potential and processes dominated by the dynam- 
ical mixing of the ISM by star formation. We at- 
tempt here to place these observations into the 
context of both disciplines, and begin by com- 
paring M101 to the well-known but not so well- 
understood ISM laboratory that is the Milky Way. 

M101 is larger than the Milky Way, having a 
£>25 about twice that of the Galaxy, but has only 
about half the mass of the Milky Way, though with 
a larger mass fraction in neutral gas (see Table 1). 
The far-infrared luminosity, a very crude measure 
of the star formation rate, is about that estimated 
for the Galaxy, implying that the star formation 
rate per disk area, or per gas mass, is substantially 
lower than the Galactic rate. Comparing M101 for 
which L fir /D 2 25 = 1.37 x 10 40 ergs s' 1 kpc 2 to 
the list of that quantity for edge-on spirals from 
Rand (1998), and the correlation with extraplanar 


diffuse ionized gas, one can see that M101 is not 
expected to have a very extended “DIG” layer. 
This is not to say that star formation in M101 is in 
any way “quiescent”; the several giant HII regions 
are more reminiscent of those in the Magellanic 
Clouds than the star forming regions seen in the 
Milky Way. 

The M101 ISM may be quite similar to that of 
the Milky Way. Kregel et al. (2002), using the 
work of Bottema (1993) has noted that the stel- 
lar velocity dispersion perpendicular to the disk is 
equivalent to the scaled circular v max . For spi- 
ral galaxies, the scale height of the old stellar 
population is nearly constant with radius (e.g., 
van der Kruit k Searle 1982). Thus, from the 
optical B-band profile (Okamura et al. 1976), the 
HI rotation curve and vertical velocity dispersion 
profile (Kamphuis 1993), and using the formulary 
of van der Kruit (1988), one can calculate the 
mid-plane density of the neutral gas, and finds 
p 0 ~ 0.35 cm' 3 for the radii of interest (4-8 kpc). 
This value is similar to the mid-plane density of 
neutral gas in the solar neighborhood (0.22-0.38 
cm" 3 , Lockman k Gehman 1991). (The calcu- 
lation uses the thin disk approximation which is 
valid for regions with a flat rotation curve; the ro- 
tation curve for M101 is gently rising at the radii, 
so the results must be viewed with some skepti- 
cism.) If one assumes that the M101 disk has a 
similar CNM/WNM/WIM structure as the Milky 
Way, then given the similar neutral gas density, 
one would expect the hot structures, if they are at 
all pressure confined, to have properties similar to 
those in the Milky Way. 

The prime difficulty with observing hot gas 
within the Milky Way is the high absorbing 
columns built up over lines of sight within the 
disk that limit our view of the Galaxy to the near- 
est .5-1.5 kpc at 0.5 < E < 1.0, and even shorter 
distances at lower energies. The filling factor of 
X-ray emitting gas within the disk is thus a hotly 
debated topic. We know of several small isolated 
examples of bubbles of ~ 10 6 K gas (Snowden 
2001), and one example of a superbubble filled 
with ~ 2 x 10 6 K gas (Egger 1993). At higher en- 
ergies, where we sample more of the Galactic disk, 
several hot Galactic ridge components have been 
identified (e.g., Kanedaet al. 1997, with kT = 0.8 
and kT = 7) but their luminosities and filling 
factors are poorly known. 
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Table 3 

Spectral Fit Parameters 


Region 

Counts* 
(0.45-2.0) keV 

Area 

arcmin 2 

N H 

10 20 cm " 2 

Ti 

kT 

EMj 
cm -6 pc 

% 

lO 20 cm - 2 

t 2 

kT 

em 2 

cm -6 pc 


V 

NGC 5461 
NGC 5462 

706/149 

1199/331 

0.607 

1.518 

oo _ 16 

.9* 

5411° 

1 ee + .038 
• IOO -.014 

261 + 034 
,zo — .219 
i 7e-f -013 
* 1/O -.028 

7 9 + 8.0 

1 '*- 5.7 

048 + 005 
.U40_ 024 

1 o+4.2 

i>S -0.9 

125 

•178 

io.o ^ 8 ° 5 

1.58 

1.23 

.968 

29 

27 

47 




.9* 

941 + 026 

no O + -006 

.Uoo_ 006 

407 1.%1 

•654^?“ 

043+036 
U4 -.023 

.938 

45 


Quantities with asterisks are fixed. All uncertainties are 90% confidence intervals. 

‘The first number is the raw counts in 0.45-2.0 keV, the second is the calculated number of particle background counts. 


About our halo we know a little more. Using 
the ROSAT All-Sky Survey, Kuntz & Snowden 
(2000) divided the Galactic halo emission into two 
components, a hard component (kT = 0.250 1'qII) 
whose intensity remains constant across the high 
latitude sky, and a patchy soft component (kT = 
0.098^ 037). (Whether these temperatures will 
survive the higher spectral resolution of the Chan- 
dra and XMM-Newton era is a topic of ongoing re- 
search.) The extraplanar X-ray emission from the 
Galactic bulge has been modeled as a cylinder of 
radius 5.6 kpc having an expo nent ial scale-height 
of 1.9 kpc and a temperature^^ = .343 (Snow- 
den et al. 1997). 

In M101, of course, we are observing both disk 
and halo gas. Assuming that half of the emission 
is behind the absorption of the galactic disk, the 
total luminosity of diffuse emission is 1.53 x 10 39 
ergs s' 1 and 0.73 xlO 39 ergs s' 1 in 0.45-2.0 keV for 
the kT = 0.2 keV and kT = 0.7 keV components. 
The total luminosity of these two components is 
^ .3% of the supernova luminosity, assuming 2 
supernovae per century (Matonick & Fesen 1997). 

5.1. Disk 

The bulk of the disk of M101 is covered by low 
surface brightness X-ray emission which traces the 
spiral arms. The low surface brightness emission 
in M101 appears to consist of two different com- 
ponents, and the constancy of the hardness ra- 
tio shows that the relative amounts of these two 
components is relatively constant for the bulk of 
the emission. This uniformity suggests that the 
sources of the two components are linked to one 


another, and to star formation. A similar two- 
component medium is seen in the spectrum of 
the actively star-forming galaxy NGC 2403 (Scd), 
where the best fit two MeKaL-component model 
has temperatures of kT = 0.18 ± 0.03 and kT = 
0.73±0.07 (Fraternali et al. 2002), while the early- 
type spiral NGC 1291 (Sa) appears to have a single 
thermal component with a temperature kT = 0.33 
keV (Irwin et al. 2002). As the diffuse emission is 
analyzed for similar galaxies in the Chandra and 
XMM-Newton archive, such two component media 
may become a mark of late-type disks. 

The softer component in M101 has a tempera- 
ture of ~ 0.2 keV. This is similar to the temper- 
ature found for superbubbles such as those in the 
Magellanic Clouds or that outlined by the Galac- 
tic Loop I. The brightest portion of the Loop 
I superbubble, a shock front seen obliquely, has 
kT = 0.17 — 0.30 and an emission measure of 0.02- 
0.04 cm -6 pc (Egger & Aschenbach 1995)). The 
emission measure of the softer component of the 
M101 disk is similar, 0.018-0.038 cm' 6 pc, sug- 
gesting a covering fraction near unity for super- 
bubbles, or, more reasonably, that a substantial 
amount of the gas is actually in the halo. An ob- 
ject similar to the Loop I superbubble (radius=150 
pc) would subtend an angle of 473 at the dis- 
tance of M101, and have a surface brightness of 
~ 3 x 10~ 3 counts cm' 2 s" 1 arcmin' 2 . Due to 
its small size, after the smoothing applied to pro- 
duce the data shown in Figure 1, it would have 
a surface brightness of only - 3 x 10' 6 counts 
cm 2 s 1 arcmin' 2 . In the presence of other dif- 
fuse emission, the Loop I superbubble would be 
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undistinguishable in the disk of M101. 

It is interesting to note that unlike NGC 2403, 
where the sound speed of the softer compo- 
nent was comfortably below the escape speed, 
in M101 the sound speed of the softer component 
(( 7 A:T//im p )5 ~ 230 km s" 1 ) is just barely com- 
patible with the escape velocity (y/2v c ~ 240 km 
s -1 ). However, this simple comparison should not 
be taken as indicating that the gas is escaping 
the galaxy. Breitschwerdt k Schmutzler (1999) 
modeled galactic chimney flows in a dynamically 
self-consistent manner and found that the adia- 
batic expansion of the gas out of the disk cools 
the gas significantly, but leaves a much higher, 
non-equilibrium ionization structure frozen into 
the gas. Their model starts with an initial CIE 
gas at 2.5 x 10 6 K, but the gas cools rapidly with 
increasing height above the disk to ~ 10 4 K, but 
the gas has a spectrum more similar to that of a 
10 6 K gas. Given that the bulk of the temperature 
information for the soft component is derived from 
the ratio of the oxygen lines, our temperature may 
not indicate the true gas temperature, but only 
the ionization temperature of the oxygen. 

Kahn (1990, 1997) has placed a more restric- 
tive constraint on this gas. Calculations for the 
Galaxy shows that the galactic chimney/fountain 
gas cools long before reaching the scale height indi- 
cated by its initial temperature, a result confirmed 
by the modeling of de Avillez (2000). Thus, the 
fountain gas is supported dynamically, does not fill 
the potential well, and does not produce a quasi- 
hydrostatic halo. 

If one does assume that the soft component has 
an exponential scale-height kT / fi ni p g , where g — 
4 ttGE, and E is the surface mass density (derived 
from the optical B-band profile under the assump- 
tion that the mass-to-light ratio is similar to that 
of the solar neighborhood, T b =2.3 MqLq ) then 
the peak electron density is 2-3x10 3 cm at 
radii of ~ 5 kpc. This value is 2 to 3 times smaller 
than the n e derived for the Local Hot Bubble 
(n e = .006 cm" 3 , kT = .108 keV), and ~ 5 times 
smaller than the densities derived either from the 
Eridion Bubble ( n e ~ .015, cm 3 , kT = .17 keV 
Guo et al. 1995) or the walls of the Loop I super- 
bubble (n e ~.01-.015 cm- 3 , Egger 1993), but is 
similar to the tenuous medium filling the center 
of the Loop I superbubble (n e = 3 x 10 3 cm 
Egger 1993). Given the softness of the Local Hot 



Fig. 13 — Dotted: The model spectrum fit to 
spectrum A-t-B. Solid: The soft and hard thermal 
(MeKaL) components. Dashed: The contribution 
by dwarf stars to the soft and hard components. 

Bubble and its overburden (typically a few x 10 20 
cm' 2 ), similar objects in M101 will not contribute 
substantially to the Chandra spectrum. Assuming 
that the soft component is formed by a population 
of objects such as the Loop I Bubble or the Erid- 
ion Bubble, and assuming that such objects have 
a signficant z extent, the filling fraction of hot gas 
would be ~ 0.2. 

If, instead, one assumes the Breitschwerdt k 
Schmutzler (1999) model, then the amount of 
truly hot gas will be a fraction of the amount in- 
dicated by the X-ray flux. Using the approxima- 
tion that the flux oc 1 /z given in Breitschwerdt k 
Schmutzler (1999) one finds about .3 of the total 
flux is due to truly hot gas. Since the fraction of 
the total flux falling within the Chandra bandpass 
may vary with z-height, a better estimate must 
wait for detailed modeling. 

The harder component has a temperature of 
~ 0.7 keV. This temperature is similar to that 
expected from the aggregate emission due to late- 
type dwarf stars {kT ~ 0.8 keV). We have used a 
model similar to that used by Kuntz k Snowden 
(2001) to calculate the emission by dwarfs (spec- 
tral types F through M) per square parsec in the 
solar neighborhood, and scaled the result by the 
disk surface mass density. It should be noted that 
while the Kuntz k Snowden (2001) model was 
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normalized to the number of K and M dwarfs in 
the solar neighborhood, and matched the number 
of Galactic stellar X-ray sources at high Galac- 
tic latitude, it did not match the mid-plane stel- 
lar mass density. Further, the Kuntz &; Snowden 
(2001) model, since it matched the number counts 
at high Galactic latitude, was not sensitive to the 
substantial uncertainties in the X-ray luminosities 
of young stars, whereas a sum through the entire 
disk is quite sensitive; even though there are one- 
tenth the M dwarfs in the 0.0-0. 1 Gyr age range as 
in the 1-10 Gyr range, their X-ray luminosities are 
ten times greater, so their contribution to the X- 
ray emission is comparable. Thus, considering the 
uncertainty with which one can apply the Galactic 
IMF and SFR to M101, the calculation of the con- 
tribution of dwarf stars to the X-ray emission in 
M101 is a rather uncertain estimate. We find that 
dwarf stars account for ^ 35% of the 0.45-2.0 keV 
band flux of the harder component, only ~ 2% of 
the softer component, and ~ 11% of the total flux. 

The temperature of the harder Chandra com- 
ponent is also similar to the u soft” Galactic com- 
ponent (kT = 0.8 keV) found by Kaneda et al. 
(1997) in their ASCA study of the Scutum arm. 
(Note that their ’’soft” component is much harder 
than the ’’hard” component of Kuntz & Snow- 
den (2000).) In their case, the spectrum was 
modeled with a non-equilibrium model requiring 
n e t = 10 9 1 cm -3 s. Replacing our hard MeKaL 
component with an NEI component does little to 
change the fit and requires n € t = 10 12 9 cm -3 s, 
that is, nearly an equilibrium plasma. From the 
fitted value of the absorbing column, Kaneda et al. 
(1997) estimate that the effective distance to the 
Galactic region producing their soft emission to be 
~ 2 kpc, and find a scale height of 100 pc. If their 
emission is produced along a 2 kpc line of sight 
then the expected emission from M101 would be 
1.9 x 10 -7 ergs cm 2 s -1 sr~‘; the observed value 
for the hard component of the A+B spectrum is 
4.7 x 10 -8 ergs cm 2 s -1 sr" 1 , only a factor of four 
smaller. If the Kaneda et al. (1997) emission is 
over a longer path, then the two values would be 
more similar. The harder component found by 
Kaneda et al. (1997) (kT ~ 7 keV) is too dim to 
be observed by Chandra at the distance of M101. 

Kaneda et al. (1997) attributed their 0.8 keV 
component to late stage supernovae expanding 
into regions of low ambient density, conjecturing 


that supernovae are abundant, one per GIS field 
of view, and reheating regions previously swept 
by older supernovae. Although the M101 equiva- 
lent is likely to be produced by the same means 
as in the Milky Way, the M101 spectra do not 
posess the evidence for non-equilibrium plasmas 
that would point to young SNR. Absence of evi- 
dence is not, however, evidence of absence, and the 
presence of such a component may be obscured 
by other emission components. A further cau- 
tion concerning the supernova attribution can be 
derived from the hydrodynamic supernova mod- 
els of Shelton (1999). In those models of indi- 
vidual supernova in low density ambient media, 
typical time-averaged temperatures were < 0.35 
keV, substantially cooler than the observed com- 
ponent. Greater ambient densities lead to signif- 
icantly higher emissivities, but not significantly 
higher temperatures. 

In M101, there are likely multiple sources for 
the hard component. Active star formation re- 
gions such as R136 in 30 Dor, for example, can pro- 
duce hot thermal components with kT ~ 0.8 keV 
(Wang 1999). A covering fraction of ,l-.2% by 30 
Dor-like objects could provide the observed emis- 
sion measure for the hard component, but would 
over-produce the soft emission by a factor of ~ 5. 

Two Component Structure: As has been men- 
tioned previously, it is not surprising that the spec- 
trum is statistically well fit by two components 
in collisional ionization equilibrium. The emission 
occupies only a small range of the energy interval 
accessible to Chandra , 0.45 to 1.0 keV, with a reso- 
lution of 0.12 keV. Within this energy interval, the 
emission is dominated by the OVII/OVIII lines at 
.560 and .650 keV, and the Fe complex at 0.7-1.0 
keV . The softer component produces the bulk of 
the oxygen flux, while the harder component pro- 
duces the Fe emission. A range of temperatures 
probably exists, but we are determining the most 
common gas temperatures. 

A similar problem was experienced by the even 
lower resolution ROSAT studies of the Galactic 
halo. Kuntz & Snowden (2000) suggested that 
the temperatures of the two Galactic halo com- 
ponents might be indicative of an isochorically 
cooling gas which passes through two regions of 
increased stability on its path down the cooling 
curve, an argument analogous to that made by 
Gehrels & Williams (1993) to explain the bimodal 
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temperature distribution of emission from RS Can 
Ven stars. 

Such an argument has a more limited appli- 
cability in this case. The cooling curve derived 
from the Raymond k Smith (1977) model al- 
lows regions of relative stability at A;X’=0.068- 
0.121 (0.79-1.4xl0 6 K) and A;T=0. 388-0.646 (4.5- 
7.5 xlO 6 K), and similar results are obtained from 
APED (Smith et al. 2001) models (Smith 2002). 
Significantly altering the metallicities does change 
the location of the inflection points, but not sig- 
nificantly for the range of metallicities observed 
in M101. The temperature of the soft component 
falls several cr above the lower inflection. The tem- 
perature of the hard component is marginally in- 
compatible with the upper limit. However, as very 
hot gas produced in the disk is not likely to cool 
isochorically, and the isobaric cooling curve does 
not have regions of enhanced stability, the tem- 
perature of the upper component is unlikely to be 
due to the cooling curve. 

5.2. Bulge 

The bulge of M101, as seen in the optical data 
of Okamura et al. (1976), can be characterized as 
having an exponential profile with a scale length of 
100 pc, much smaller than the bulge of the Galaxy, 
which has a scale length of ^ 1 kpc. In M101 the 
extended X-ray nucleus has a scale length similar 
to that of the optical bulge, and the bulge-to-disk 
ratio in the X-ray is similar to that of the optical. 
Given the small size of the bulge and the scaling 
of the X-ray emission with optical light, it is per- 
haps not surprising that, unlike the Galaxy, the 
X-ray emission does not show a marked change in 
temperature between disk and bulge. In the X- 
ray range, there are no signs of nuclear activity, 
and the nuclear point sources are consistent with 
ordinary X-ray binaries. 

5.3. Extended Halos 

Since the classic paper by Spitzer (1956), the 
high energy model of galaxies has included a 
nearly spherical hot corona or halo, either in hy- 
drostatic equilibrium, or in some similar state. 
The existance of such a halo seems reasonable 
given that the virial temperature of a galaxy like 
the Milky Way is ~ mv\/ ( lk = 2 x 10 6 K, where 
v c is the circular velocity. Models of the soft 



Fig. 14. — The radial profile for position angles 
60° to 110°. The solid line excludes the giant HII 
regions; the dashed line includes the giant HII re- 
gions. 

X-ray emission from our own Galaxy usually in- 
clude such a component (Wang 1998; Pietz et al. 
1998), though Kuntz k Snowden (2001) suggested 
that the rather featureless component that domi- 
nates the ROSAT § keV emission may not be due 
to a Galactic corona, but rather to cosmological 
sources. 

Observations of other galaxies in the ROSAT 
era were generally equivocal; the only extended 
halos detected around non-starburst, edge-on 
galaxies were arguably due to high star forma- 
tion rates (NGC 4631 and NGC 891). Benson 
et al. (2000), noting that simple cooling flow mod- 
els of galaxies predicted Lx oc searched for 
extended X-ray halos around three massive galax- 
ies. Although these galaxies were expected to 
have L x ~ 10 42 ergs cm' 2 s“\ the observed 95% 
upper limits were 3.8, 1.2, and 0.4 xlO 41 . 

Toft et al. (2001), executing SPH models of 
galaxy formation from z ~ 20 — 40 to the present, 
found their galaxy halos to be at least an order 
of magnitude fainter than the halos produced by 
the simple cooling flow models of galaxy forma- 
tion. M101 is voluminous, but not massive. The 
L x — v c relation found by Toft et al. (2001) sug- 
gests that M101 {v c ~ 170 km s' 1 ) should have a 
halo with Lx = 1 “ 3 x 10 38 ergs s" 1 in the 0.2- 
2.0 keV band, T < 10 6 AT, and the X-ray emission 
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should be fairly closely confined to the disk. 

Figure 14 shows the radial profile of M101 be- 
tween position angles 60° and 110° (measured 
from North). At radii greater than ~ 7 kpc, the 
observed surface brightness drops sharply to the 
value of the background measured by ROSAT ; the 
2<r upper limit is 9 x 10 -16 erg cm -2 s -1 in 0.45- 
1.0 keV. Thus any axisymmetric halo must make 
an insignificant contribution beyond this radius. 
The Toft et al. (2001) value for halo luminosity 
is 7-20% of the measured 0. 2-2.0 keV luminosity 
of the diffuse emission within 7'. Since the X- 
ray luminosity follows the optical luminosity (i.e., 
Figure 4), and there is no excess X-ray emission 
within 7 other than the bulge, the true strength 
of an X-ray halo not due to disk processes is likely 
to be much smaller than the 7-20% predicted by 
Toft et al. (2001). From the viewpoint of ISM 
studies, this result is the anticipated consequence 
of galactic fountains. 

5.4. Flotsam 

Nucleosynthetic studies of the Milky Way sug- 
gest that galaxies continue to accrete material, 
~ 1M 0 year -1 for the Galaxy (Alibes et al. 2001; 
Timmes et al. 1995), and generally cite the (as 
yet unknown) mass of high velocity clouds as a 
plausibility argument. The luminosity due to the 
infalling material should be ~ 10* ergs s ~ 1 , five 
times the observed amount of diffuse X-ray emis- 
sion in this Chandra observation. This point may 
not be consequential, as estimates for the current 
infall rate required for nucleosynthesis appears to 
have at least a factor of three uncertainty (e.g. 
compare Alibes et al. 2001 with Timmes et al. 
1995). Further, given the structure of high veloc- 
ity clouds, the infall is likely episodic. However, 
M101 appears to be abundant in high velocity gas, 
having ~ 2 x 10 s M 0 (Kamphuis 1993). Unfortu- 
nately, this gas is situated entirely outside of the 
field of view of this Chandra observation; XMM- 
Newton may be able to determine to what low level 
the high velocity cloud/halo interface shock con- 
tributes to the diffuse X-ray emission. 

6. Summary 

The bulk of the truly diffuse emission from 
M101 traces the spiral arms and is correlated 
with Ha and FUV emission, implying that it is 


due to ongoing star formation: hot winds from 
massive young stars, supernovae, and superbub- 
bles. The spectrum of the diffuse disk emis- 
sion is characterized by two thermal components 
with temperatures kT = 0.20 keV and k.T = 
0.75 keV. The ratio of the emission measures is 
roughly constant across the face of the galaxy with 
EM so f t /EMh a rd ~ 3. Comparison of the tem- 
peratures and emission measures of these com- 
ponents to diffuse features in the Galaxy or the 
Magellanic clouds suggest that a significant frac- 
tion of the gas forming the soft component is ex- 
traplanar. The hard component in M101 is sim- 
ilar to the soft Galactic ridge component seen in 
the Milky Way, but has a substantial contribution 
from dwarf stars. 

We do not find evidence for an axisymmetric 
halo tracing the dark matter potential. At a ra- 
dius of ~ 7 kpc, the deep interarm region shows 
no sign of X-ray emission above the level of the 
background. Given that M101 is a fairly low mass 
galaxy, the lack of such a halo may not be surpris- 
ing. 

Whether M101 has a true nuclear source re- 
mains unknown. The diffuse emission from the 
highest surface brightness portion of the nu- 
cleus has a spectrum consistent with thermal 
Bremsstrahlung absorbed by about half of the 
total M101 column density in that direction. The 
remainder of the nuclear emission has a spectrum 
consistent with that of the diffuse disk emission. 

The giant Hll regions NGC 5461 and NGC 5462 
are composed of truly diffuse emission whose spec- 
tral shape is similar to that of the disk emission, 
but the total number of counts is too small to allow 
detailed spectral fitting. The non-thermal compo- 
nent in these spectra appears to be small. 
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